The three-dimensional, three state Potts model in a negative external field 
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We investigate the critical behaviour of the three-dimensional, three state Potts model in presence 
of a negative external field h, i.e. disfavouring one of the three states. A genuine phase transition is 
present for all values of \h\, corresponding to the spontaneous breaking of a residual Z2 symmetry. 
The transition is first /second order respectively for small/large values of \h\, with a tricritical field 
fttric separating the two regimes. We provide, using different and consistent approaches, a precise 
determination of hmc', we also compare with previous studies and discuss the relevance of our 
investigation to analogous studies of the QCD phase diagram in presence of an imaginary chemical 
potential. 
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I. INTRODUCTION 

Potts models [H have been often considered in litera- 
ture to reproduce the critical properties of more complex 
physical systems 043- In the present work we are in- 
terested in the 3-state Potts model defined on a three di- 
mensional cubic lattice. The generic g-state Potts model 
is defined by the following partition function 



HAI) 



(1) 



where the spin variable <7j lives on lattice site i and can 
take q different possible values, e.g. <7j £ {0, 1, ...,(<? — 
1)}, while f3 = l/(fcsT). E and M denote respectively 
the energy and magnetization with respect to a chosen 
reference spin value a (e.g. a = 0): 



M = 



(2) 
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where J is the coupling constant, H is an external applied 
field and in the sum denotes all pairs of nearest 

neighbor lattice sites. In the following, as usual, we shall 
set J — 1 and make use of the normalized magnetic field 
h = pH. 

For h — the system has an exact symmetry, corre- 
sponding to all possible global permutations of the q spin 
values, i.e. the symmetry group is the group of permuta- 
tions S q . Such symmetry gets spontaneously broken be- 
low a given critical temperature, where the spin variables 
align themselves along a given direction. The correspond- 
ing phase transition is first order, in three dimensions, for 
q > 3, and second order for q = 2 (the system coincides 
with the Ising model in this case). For q — 3 the critical 
temperature is given by p c (h — 0) — 0.550565(10) [8(. 

The critical properties of 3D Potts models with q = 2 
or 3 at h = 0, have often been associated with those met 



at the finite T phase transition of QCD (with 2 or 3 col- 
ors) in the pure gauge limit, via the well known Svetitsky- 
Yaffe conjecture Q. The symmetry group which is spon- 
taneously broken in the high T, deconfined phase of 
SU(N) pure gauge theories is that associated with center 
symmetry, Zn, corresponding to local gauge transforma- 
tions which are periodic in the Euclidean time direction 
only up to a global group element belonging to the center 
of the gauge group. The symmetry group coincides with 
the permutation group for N = 2, while for N = 3 one 
has to add charge conjugation to center transformations 
to recover the full permutation group 5*3. The corre- 
sponding order parameter in pure gauge theories, playing 
the role of magnetization and signalling the spontaneous 
breaking of center symmetry in the deconfined phase, is 
the Polyakov loop, i.e. a closed parallel transport in the 
Euclidean time direction. 

For ft, 7^ the symmetry S q is explicitly broken to 
a residual S q -\, corresponding to permutations among 
spin values other than the chosen direction a. The case 
h > 0, in which alignment of spin variables along a is 
favoured, has been extensively studied in the literature: 
one can still distinguish two phases in which the system 
is more (low T) or less (high T) aligned along a, however 
in both phases the residual S q -i symmetry stays unbro- 
ken, so that no real phase transition is expected a priori. 
However, for q > 3, the first order transition present 
at h = persists also for non zero positive values of h, 
till a critical endpoint is met, after which the transition 
disappears. Such critical endpoint is expected to be in 
the Ising 3D universality class and for q = 3 it has been 
located at (p c ,h c ) = (0.54938(2), 0.000775(10)) HOJ]. 

In the analogy with finite T pure gauge theories, the 
case h > corresponds to adding dynamical fermions 
of mass to and in the fundamental representation of the 
gauge group (the limit h — >• corresponding to m — > 00): 
that induces an effective coupling to the Polyakov loop, 
which breaks center symmetry and aligns the Polyakov 
loop along the positive real direction, while the residual 
charge conjugation symmetry stays unbroken for all val- 
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ues of T. 

The case h < is quite different. Indeed in this case 
the reference state a is disfavoured, so that low T order- 
ing happens along one of the remaining q — l states: there 
is therefore an order /disorder transition associated with 
the spontaneous breaking of the unbroken S g -i symme- 
try group. Actually, in the limit of large \h\, the system 
becomes completely equivalent to a (q — l)-state Potts 
model at zero magnetic field, since the disfavoured state 
disappears from the statistical ensemble. Therefore a 
true phase transition is expected for every value of \h\, 
coinciding with the transition of q- or (q — l)-state Potts 
model in the limit of zero or infinite field respectively. 

In the present paper we shall discuss the case q = 3 
in three dimensions, which is particularly interesting (as 
well as the case q = 5 in two dimensions), since in this 
case the transition at h = is first order, while the tran- 
sition at h = — oo is second order in the 3D Ising uni- 
versality class. Hence the expectation is that the first 
order continues for small values of \h\, until a tricritical 
point /itric is met, governed by mean field indexes, af- 
ter which the transition becomes second order in the 3D 
Ising universality class. An accurate verification of this 
scenario and the precise location of the tricritical point is 
the aim of our study. Notice that, for small \h\, we expect 
an interesting example of system which may be naively 
believed in the Ising 3D universality class because of sym- 
metry reasons (the relevant symmetry being Z2), but has 
instead a first order transition because of the interplay 
with different dynamical degrees of freedom, correspond- 
ing to the disfavoured state a in this case. 

Going back to the correspondence with the critical 
properties of SU (3) lattice gauge theories, switching the 
sign of h is like turning the boundary conditions of dy- 
namical fermions in the temporal direction from anti- 
periodic to periodic: fermions are not thermal any more 
and the Euclidean temporal direction can then be viewed 
as a compactified spatial direction; the effective coupling 
to the Polyakov line changes sign, so that the Polyakov 
line tends to align along one of the complex center ele- 
ments below a given compactification radius, thus break- 
ing spontaneously the residual charge symmetry (see e.g. 
Refs. [T7l - [20| ] for early lattice studies of the associated 
transition, which has been studied in the context of ori- 
entfold planar equivalence [2l], HI])- Alternatively, one 
can interpret the system as the usual thermal theory in 
presence of a purely imaginary quark number chemical 
potential such that Im(/z)/T = 7r: in that case the Z2 
breaking transition is interpreted as the endpoint of the 
high T Roberge- Weiss (RW) transitions which are met in 
the r-Im(fi) plane [23[. 

The importance of this Z2 transition and of its or- 
der for the general features of the QCD phase diagram 
has been discussed extensively in recent literature (23 - 
l29l ]. In particular, its order has been investigated by 
lattice simulations in QCD with two degenerate flavors 
in Ref. 12611 . and more recently also for the three flavor 
theory [27]: in both cases one finds a non-trivial phase 
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TABLE I: Critical exponents (see e.g. [33, [32 



structure, with the transition being first order both for 
small and high quark masses, and second order in the 
middle. Such phase structure can be mapped to that of 
the Potts model with a negative magnetic field, which is 
the subject of our study, on the large mass side; on the 
other hand, on the small quark mass side, chiral degrees 
of freedom come into play, requiring a different effective 
model description. In the context of the investigation of 
the QCD phase diagram, it is of course particularly im- 
portant to give precise estimates of the tricritical values 
of the quark mass, separating the second order from the 
first order regions. 

The study of the 3D three-state Potts model in a neg- 
ative magnetic field can be placed in the more general 
context of studies of the same model in complex mag- 
netic fields [ll|, EH H3 1 aimed at mimicking the dynam- 
ics of QCD in presence of a quark chemical potential, 
which have also considered the properties of the tricritical 
point [27| . Our purpose is that of performing a detailed 
study of the critical behavior of the system as a function 
of h, with the specific aim of determining the location 
of the tricritical field h tr i c . We will make use of differ- 
ent and consistent approaches in order to do that: the 
strategy developed for this model and the corresponding 
results can then be taken as a guideline for the analogous 
determination of the tricritical masses for the endpoint 
of the Roberge- Weiss transition in QCD [3(| • 

The paper is organized as follows: in Sec. [II] we present 
and discuss the different strategies used to investigate the 
critical properties of the system; in Sec. Mil we present 
our numerical results and finally, in Sec IIV1 we give our 
conclusions. 



II. OBSERVABLES AND NUMERICAL 
ANALYSIS SETUP 

Natural observables for the Potts model are the energy 
E, which is defined in Eq. @, and the magnetization. As 
for the latter, we replace it by a new quantity P which, 
in the analogy with QCD, plays the role of the average 
Polyakov line. In order to define P, let us associate with 
each spin variable a complex number on the unit circle 
as follows 



exp 



then we define 



p=-Y,< 



(4) 



(5) 



where V = L z is the lattice volume. Assuming that the 
state coupled to the magnetic field is a = 0, the residual 
Zi symmetry of the model corresponds to an exchange 
of the states 1 and 2, i.e. to complex conjugation for the 
complex spin variables Si and for P. Therefore, while 
E and Re(P) are even under the residual Z2 symmetry, 
Im(P) is odd and plays the role of the order parameter 
for the realization of this symmetry. 

The purpose of our investigation is that of determining 
the location of the phase transition, its order and univer- 
sality class, as a function of the magnetic field h, which 
is taken to be negative. In order to do that, we shall 
consider at first the susceptibility of the order parameter 

X^L 3 ((Im(P) 2 )-{|Im(P)|) 2 ) (6) 

and the specific heat of the system 

C = L 3 ((E 2 ) - (E) 2 ) . (7) 

The scaling of the two quantities around the phase tran- 
sition, as a function of the size of the system, is fixed by 
the respective critical indexes 

X = m v MtL 1/v ) (8) 

and 

C = Co + L a ' v h(tL 1/u ) , (9) 

where Co is a regular contribution and t = (T — T c ) /T c 
is the reduced temperature. In the following we will be 
interested in particular in the scaling of the- height of the 
peaks, which is regulated by j/u and ajv respectively, 
and in the scaling of the width of the peaks, which is reg- 
ulated by 1/v in both cases. The critical indexes which 
are relevant to the different possibilities which may take 
place {i.e. first order, second order in the 3D Ising univer- 
sality class and mean field tricritical) are listed in Tablefl] 
Another interesting quantity is the modulus of P, how- 
ever it takes contribution both from the order parameter 
and from the spin state coupled to the magnetic field, 
which is Z2 even, hence its susceptibility is expected to 
be the mixing of different contributions scaling like \ or 
C respectively, apart from the limit h — > 00, in which case 
the contribution of the state coupled to h is completely 
suppressed and the system can be mapped exactly to a 
3D Ising model. We shall not consider such quantity in 
the following. 

It is interesting to notice that, while first order scal- 
ing and 3D Ising scaling are expected to take place for 
a continuous range of values of h, tricritical scaling is in 
principle expected only for a specific value /itric, which 
we want to determine, at the boundary between the first 
order and the second order region. However the cor- 
rect expectation is to have tricritical scaling regulating a 
neighborhood of /itric, with the size of the neighborhood 
going to zero as L —> 00. Putting the question the other 
way around, we expect to need increasingly large volumes 
to discriminate between first order and 3D Ising second 
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FIG. 1: On finite volumes tricritical scaling is expected to 
dominate a finite range of h values around /itric, which shrinks 
to zero as L -> 00. 

order as we approach the tricritical field /i tr ic from cither 
side, since a fictitious tricritical scaling will mask the true 
thermodynamical limit for not large enough volumes (see 
Fig. [T] for a graphical representation of that). 

This reasoning can be made more quantitative by use 
of the so-called crossover exponents: in the thermody- 
namical limit, the true critical behaviour of the system 
can be seen only for \t\ < p 1 ^, where t is the reduced 
temperature, p is the parameter that controls the change 
of critical behaviour and <j> is the crossover exponent (see 
e.g. (32H341 ] ) . which is by definition <f> — y p /yt (yt and 
y p are the renormalization group eigenvalues of the rel- 
evant variables t and p). In our case p oc h — fttric an d 
0=1/2 [35[. On a finite lattice of typical size L, \t\ can 
be traded for L~ x l v and the previous condition becomes 
L ^ \h — htric\~ v ^] in particular, according to the known 
tricritical indexes in Table HI one expects tricritical be- 
haviour to dominate up to a critical size 

L c ~A\h — /itricT 1 (10) 

where A is some unknown constant which may be dif- 
ferent on the first order and on the second order side; 
a numerical check of this behaviour will be reported in 
Sec. MI Al (see in particular Fig. [T0|) . That implies that 
a correct and precise determination of /itric may be quite 
difficult if one looks at the finite size scaling of suscepti- 
bilities or other quantities alone. 

As an alternative and easiest way to determine /itric j 
we shall determine quantities which give a measure of the 
strength of the first order transition, such as the latent 
heat or the gap of the order parameter at the transition, 
and study the variation of such quantities as a function 
of h, in order to extrapolate the point /itric where they 
vanish, i.e. where the first order disappears, without the 
need of making simulations very close to /itric- 

We shall consider in particular the Binder-Challa- 
Landau cumulant (36| of the energy, which is defined as 
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FIG. 2: Probability distribution of the energy density for dif- 
ferent lattice sizes L and h — —0.0025, where the transition 
is first order. 




FIG. 3: Probability distribution of the energy density for dif- 
ferent lattice sizes and h = —0.01, where the transition is 
second order. 



B 4 = 1 - (E 4 )/(3(E 2 ) 2 ). It can be shown (see e.g. [33]) 
that near a transition B4 develops minima whose depth 
scales as 
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0(L- 3 ) 
+ 0(L- 3 ) (11) 
- E- and e = 



E-). In particular the thermodynamical limit 
is less than 2/3 if and only if a latent heat is 
present; to simplify the notation in the following we will 
use the shorthand B = 3 — Bi\ min . 

A different, but analogous quantity is the gap of the 
order parameter, A, which can be extracted by looking 
at the scaling of the maximum of its susceptibility, x> 
and using the relation, valid in the large volume limit for 



a first order transition, 



Xn 



T 3 

A+^-A 2 . 
4 



(12) 



Both A^ and A are expected to vanish as we approach 
the tricritical field /itric from the first order side. In par- 
ticular, the leading order expected behaviour is the fol- 
lowing (see or [38[ for a brief summary) 



and 



A E oc \Jh- h t 



A oc y/\(h- h t nc)\og(h - /l t ric)| 



(13) 



(14) 



Another useful quantity is the fourth-order cumulant 
of the order parameter. This is usually defined by 
(M 4 )/(M 2 ) 2 ([H, S3), where M is the order parameter, 
and is typically used in the study of second order transi- 
tions. Since in this work we will analyze mainly region of 
the parameter space in which first order transitions are 
present, the connected form 



U A = 



mm 

({SM) 2 ) 2 



5M = \lmP\ - (\ImP\) (15) 



appears to be best suited to disentangle the fluctuations 
inside a thermodynamical phase from the tunneling be- 
tween the two sectors with different Z2 magnetization. 
Another reason to prefer the connected form is that at 
the critical point it develops a minumum, making thus 
possible to obtain the value of the cumulant at the tran- 
sition, without introducing cross-correlations with other 
observables, or between different lattice size data. As a 
last point, we note that the relative error of the cumulant 
value at transition turned out to be smaller by a factor 
2 for the connected cumulant than for the usual one. 

For a first order transition it is simple to show that 
in the thermodynamic limit U4 — > 1, by using a double 
gaussian approximation for the distribution of the order 
parameter. For a second order transition it can be shown 
that the value of U4 at the transition is a renormalization 
group invariant ( [3QL [ioj ) , so that the intersection point 
of U4 calculated on two lattices of different size can be 
used as an estimator of the transition point. It can also 
be shown that the slope of U4 at the transition point V\ 
satisfies the relation 



dU 4 (bL) 



dUi{L) 



= b 1 / v 



(16) 



thus giving an estimate of the v critical index. 



III. NUMERICAL RESULTS 

The first order transition, which is already quite weak 
at h = 0, gets weaker for negative h values, so that we do 
not need to use algorithms specifically designed for strong 
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first orders, like the multicanonical one. While approach- 
ing the tricritical point autocorrelation times grow up, 
however, since we will perform our simulations mainly in 
the first order region, this slowing down is not expected to 
be too significant for our study. Numerical simulations 
have thus been performed using a standard Metropolis 
algorithm. 

Collected statistics have been of the order of 10 7 -j- 10 s 
independent configurations for all volumes and parame- 
ter sets explored; numerical simulations have been per- 
formed on GRID resources provided by INFN. 



A. Discerning the critical behaviour from finite 
size scaling. 

One way to discern between a first order and a second 
order critical behaviour 1 is to look at the distribution 
of physical observables, like the energy, at the transi- 
tion point: that is expected to develop a double peak 
structure, in the thermodynamical limit, for a first or- 
der transition, while it stays single peaked in the second 
order case. In Figs. [2] and [3] we show two examples, for 
h = —0.0025 and h = —0.01 respectively, where the sit- 
uation is quite clear: h = —0.0025 clearly belongs to the 
first order region, while h = —0.01 appears to be on the 
second order side. 




L 

FIG. 4: Scaling of the specific heat peak with L for h — 
-0.0025. x 2 / d -°- f - - 0.4 (range of fit: L > 50). 

Such conclusions are confirmed by looking at the scal- 
ing of the height of the specific heat peak. For h = 
—0.025 (see Fig. @| a cubic term in L, which is char- 
acteristic of first order, nicely fits the behaviour on the 
larger volumes. For h = —0.01 (see Fig. [5]) the situation 



For the sake of simplicity we will speak of "critical behaviour" 
also for the case of first order transitions, although this is not 
completely appropriate. 
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L 

FIG. 5: Scaling of the specific heat peak with L for h = —0.01. 
In this case the correct scaling with 3D Ising critical indexes 
is visible already from moderate size lattices. xVd.o.f. ~ 0.95 
(range of fit: L > 20) 
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FIG. 6: Scaling of the specific heat peak with L for h = 
—0.005. According to our determination of Titric, the system 
belongs to the 3D Ising universality class, however data scale 
according to tricritical indexes, with small deviations appear- 
ing only on the largest available volumes, x 2 /d.o.f. ~ 1.4 
(lattices with L > 100 are not included in the fit). 



is also quite clear and data correctly scale according to 
3D Ising critical indexes. 

What is less clear is the critical behaviour of the spe- 
cific heat peak for h = —0.005, which is shown in Fig. [5] 
Data scale linearly with L, i.e. according to tricritical 
indexes, for a large range of lattice sizes, with small de- 
viations, going in the direction of a smallest value of a/v 
(hence in the direction of the 3D Ising class), visible only 
on the largest sizes explored, L > 100. Our subsequent 
analysis will clearly show that for this value of h the sys- 
tem belongs to the 3D Ising universality class, however 
it would have been difficult to state that clearly from the 
scaling of the specific heat alone: tricritical indexes reg- 
ulate the system behavior till L ~ 100, completely mask- 
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ing the correct thermodynamical limit, which would be 
evident only on much larger lattices. We expect the sit- 
uation to be worse and worse as one gets closer to /itric 
(see the previous discussion in Sec. ITT|) . 
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FIG. 7: Scaling with L of the peak of the order parame- 
ter susceptibility for h — —0.0025. The contribution of two 
different terms, corresponding respectively to first order scal- 
ing and tricritical scaling, is needed to correctly fit our data. 
x7 d -°-b ~ 0.45 (range of fit: L > 50). 




FIG. 8: Scaling with L of the peak of the order parameter 
susceptibility for h = —0.005. In this case it is hard to disen- 
tangle tricritical from 3D Ising scaling, since j/v = 2 in the 
first case and y/v ~ 1.968 in the second case. xVd.o.f. ~ 1.4 
(range of fit: L < 130). 

The situation is even more difficult when studying the 
scaling of the peak of the order parameter susceptibility, 
X, since in this case the relevant critical index, j/v, prac- 
tically coincides for the tricritical ("f/f = 2) and 3D Ising 
cases ("f/u ~ 1.963). 

For h — —0.0025 one clearly sees a first order contri- 
bution (see Fig. [7]). Notice however that, in order to cor- 
rectly fit data, it is necessary to take into account also a 
small but non-zero contribution proportional to L 2 ; this 
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FIG. 9: Scaling with L of the half height width of the peak of 
the order parameter susceptibility for h = —0.005. While till 
moderate size the scaling is compatible with a tricritical one, 
on larger volumes a deviation from the tricritical behaviour is 
clearly seen. xVd.o.f. ~ f.3 for the small volumes (L < 100), 
X 2 /d.o.f. ~ 0.3 for the larger ones (L > 80). 



is the dominant term in the case of tricritical scaling, 
therefore we can interpret that as evidence for a non- 
negligible influence from a possibly close tricritical point . 
The fit with the functional form a + bL 3 + cL 2 gives the 
estimates for the parameters a — 710(70), b — 0.0139(3) 
and c = —0.37(3), meaning that the tricritical correc- 
tions to first order is about 20% on the largest lattices 
explored for this h value. For h = —0.005 there is no 
hope to discern between mean field and 3D Ising (see 
Fig.©. 

A better probe in this case is furnished by the width of 
the susceptibility peak at half height, which is expected 
to scale like L -1 / 1 ', since \ jv changes appreciably from 
mean field tricritical to 3D Ising (see Table H]). Data for 
the width are shown in Fig. |H1 two different regimes are 
visible, the first for L < 100 regulated by the tricritical 
exponent (y = 1/2) and the second, for larger lattices, 
where this index is sensibly larger, v = 0.526(3). We 
explicitly note that this value for v is to be regarded 
just as an "effective", size dependent, index interpolating 
between the tricritical {v = 1/2) and Ising one (y = 
0.63). That confirms what already found by looking at 
the scaling of the specific heat. 

Let us summarize and comment the results contained 
in this subsection. Discerning the correct critical be- 
haviour from the finite size scaling analysis of susceptibil- 
ities or other quantities may be a difficult task since, as 
expected, tricritical behavior masks the correct asymp- 
totic scaling behaviour for some range of lattice sizes, 
which increases as we get closer to the tricritical field 
/itric according to the tricritical crossover exponents, as 
summarized in Eq. (|10[) . We have tried to verify quantita- 
tively the prediction reported in Eq. (|10p by estimating, 
for each value of h on the first order side, the critical size 
L c such that for L > L c the scaling of the maxima of 
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the specific heat is well described by a first order scal- 
ing. Results are reported in Fig. [TUJ Eq. (fTU|) is well 
verified by using the value of fttric obtained and reported 
in Section IfflTfl 



-J 80 



FIG. 10: On the vertical axis it is plotted the size L c such 
that for L > L c the scaling of the maxima of the energy 
susceptibility is well described by a first-order scaling. The 
red line is a fit of the form a + b/\h — fttric| where frtric = 
-0.00415(3) (see Section HUB) . 



The difficulties are generally larger on the second or- 
der side than on the first order one, and we can easily 
understand why: the growth of susceptibilities is larger 
for mean field tricritical indexes than for 3D Ising crit- 
ical indexes, hence a fake tricritical behavior can mask 
3D Ising indexes for a large range of lattice sizes; on the 
other hand a first order behaviour, which implies a faster 
growth of susceptibilities with respect to the tricritical 
one, is in general more easily detectable as a correction 
to tricritical behaviour starting from smaller lattice sizes. 
As a last comment, we note that such difficulties make 
it preferable to look at the scaling of the specific heat 
rather than at that of the order parameter, since the 
critical index regulating the growth of the specific heat 
with L, a/u, changes more drastically when going from 
first order {a/v = 3) to tricritical {a/v — 2) and to 3D 
Ising (av ~ 0.175). 



B. The latent heat, the order parameter gap and 
determination of /i tr ic 

As explained in Section [TTl we will now determine the 
parameters which fix the strength of the first order tran- 
sition taking place for small values of \h\, in order to 
extrapolate the critical value /itric at which the first or- 
der disappears. The parameters are the latent heat, or 
equivalently the minimum of the Challa-Landau-Binder 
cumulant defined in Eq. (fTT|) . and the gap of the order 
parameter, which can be extracted from the large volume 
limit of the maximum of its susceptibility x, see Eq. (|12p. 

In Fig.QTJwe show the quantity B (see Eq. (|TT|) and the 



related discussion) as a function of 1/V for h = —0.0025. 
It clearly extrapolates to a non-zero value for V — > oo, 
a = 8.35(4) x 10~ 4 , with both 1/V and 1/V 2 corrections 
visible in the range of explored volumes. For the same 
value of h and using the same fit shown in Fig. from 
the coefficient of the cubic term in L we extract the order 
parameter gap, A 2 = 1.689(5) x 10 -3 . In Fig. [T2"l instead 
we show the case h = —0.005, together with a power law 
fit B oc L a . If we try instead B = bo + L a we get for &o the 
result consistent with zero shown in Tab. HH indicating 
that no latent heat is present. 




FIG. 11: Binder-Challa-Landau cumulant of the energy for 
h = —0.0025 > /itric- B extrapolates to a non zero value as 



V 



oo. x /d.o.f. 



1.4. 




0.01 0.02 0.03 
1/L 



130 120 110 100 90 80 70 60 50 40 

1/L 

FIG. 12: As in Fig. [TTl for h = -0.005 < /itric B goes 
to zero but the power law changes beyond a given size sep- 
arating tricritical from 3D Ising scaling: a = —1.77(1) and 
b = -1.98(1). x7 d -°-f- ^ 1-4 (L < 100) and x 2 / d -°- f - - 0.4 
(L > 80) respectively. 

We have applied the same procedure to all values of h 
where the first order transition is clearly detectable on 
the explored volumes, obtaining the values for B and the 
gap reported in Table ITT] From those values, and using 
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h 


B 


A 2 


-0.002 


1.17(2) x 10~ 3 


1.45(2) x 10~ 2 


-0.0025 


8.29(6) x 10~ 4 


1.28(7) x 10~ 2 


-0.003 


5.1(1) x 10~ 4 


7.8(3) x 10~ 3 


-0.0035 


2.78(1) x 10~ 4 


5.6(3) x 10~ 3 


-0.0038 


1.7(2) x 10~ 4 


2.7(4) x 10~ 3 


-0.005 


5(11) x 10~ 6 


-2(2.8) x 10" 4 



TABLE II: Estimated values for the thermodynamical limit 
of B and A 2 . 



the expected behaviors reported in Eqs. (1131) and (T141) . 
we can fit the value of /itric from both quantities. Results 
are reported in Fig. [T3] we obtain htiic = —0.00410(5) 
from the extrapolated minimum of the cumulant, and 
^tric = —0.00412(7) from the order parameter gap. The 
two values are in perfect agreement with each other and 
with the outcome of the finite size scaling analysis re- 
ported above; however the finite size scaling analysis 
alone would have not been able to locate /itric with such 
precision. We also notice that our determination for 
h trlc is in good agreement with the results reported in 
Refs. [H[ and [13] (Fig. 5 in both references), whose es- 
timate was 2 fttric = -0.00445(20). 




FIG. 13: Extrapolation of the Binder cumulant and of the 
order parameter gap in order to extract the tricritical value 
of the magnetic field /itrio 

As an alternative, independent way to locate h tr - lc , we 
have studied the cumulant U4 defined in Eq. (fT5")) . The 
theoretical expectation is that increasing the lattice size 
U4 — > 1 for h > fttric, U4 y Uf lc for h = /i tr i C and 
U4, — > U^ ins for h < /itric- In particular the cumulants 
calculated on different lattices arc expected to intersect 
at the tricritical point, with a slope increasing as L 1 ^ 
(see Eq. [HI), with v = 1/2. 



Numerical results for U4 are reported in Fig. [14] and 
Table IIIII the location of the intersection poin t is de- 
termined by using the method exposed in [42|], §111. B: 
a scaling law of the form U4 = f((h — ht r i c )L v ) is as- 
sumed and, since we are sufficiently close to the tricrit- 
ical point, we can develop f[x) in power series around 
x = (also scaling corrections are usually to be taken 
into account, see the discussion in (42j); a fit is then per- 
formed in the expansion parameters taking into account 
the data measured at different L values. By using data 
for -0.005 < h < -0.002 and L > 40 the fit has 19 
d.o.f. and x 2 /d.o.f. ~ 0.8 The estimated location of the 
tricritical point is h tr i c = —0.00415(3). In Fig. [15] it is 
shown that the derivative of the cumulant scales with the 
expected critical index. 




-0.005 



h 



FIG. 14: Binder cumulant of the order parameter at the criti- 
cal temperature, as defined in Eg. 1151 The curves at different 
volumes intersect at /itric- 



-200 1 L- 

40 2 



50 2 



«r 



7() : 



sir 



FIG. 15: Slope of the Binder cumulant of the order parameter 
at the critical temperature. The line is a linear fit. 



2 P. de Forcrand, private communication. 



h 


L = 40 


L — 50 


L — 60 


L — 70 


L = 80 


-0.0015 


1.496(2) 


1.410(2) 


1.325(2) 


1.244(2) 




-0.002 


1.573(4) 


1.512(5) 


1.435(4) 


1.367(2) 


1.301(3) 


-0.0025 


1.648(4) 


1.598(4) 


1.542(5) 


1.496(4) 


1.431(5) 


—0.003 


1.710(2) 


1.677(2) 


1.646(3) 


1.612(3) 


1.573(4) 


-0.0035 


1.770(3) 


1.745(4) 


1.739(7) 


1.719(6) 


1.691(7) 


-0.0038 


1.794(4) 


1.781(5) 


1.782(6) 


1.769(5) 


1.761(6) 


-0.005 


1.880(4) 


1.906(4) 


1.913(7) 


1.932(6) 


1.955(5) 


-0.01 


2.080(6) 


2.071(6) 


2.121(7) 


2.116(8) 


2.123(9) 



TABLE III: U4 values at the transition for different lattice 
sizes and magnetic field. 



h 




-0.0015 


0.549537(4) 


-0.002 


0.549237(2) 


-0.0025 


0.548942(1) 


-0.003 


0.548652(1) 


-0.0035 


0.548358(3) 


-0.0038 


0.548199(2) 


-0.005 


0.5475152(6) 


-0.01 


0.545071(8) 


-0.5 


0.484166(5) 


-1.0 


0.465188(4) 


-1.5 


0.45576(1) 


-2.0 


0.450591(4) 



TABLE IV: Estimated values for j3 c at fixed h. 



C. Critical temperature 

We conclude the presentation of our numerical results 
by analyzing the behaviour of the critical temperature 
as a function of h, j3 c (h). Our determinations of /3 C are 
summarized in Table llVl and they have been obtained by 
using the number of phase criterion for first order tran- 
sitions ([§[) , while for the second order ones the crossing 
point of the order parameter cumulant was used (see e.g. 
40]). 

We expect that for large negative values of h the state 
coupled to the magnetic field disappears from the system 
dynamics, which then becomes completely equivalent to 
that of a 3D Ising system. That must be visible from the 
behaviour of /3 c (h) which should approach two times 3 
/3 C (Ising) = 0.2216546(10) (0) as h -> -00. In Fig. [IB] 
we show the quantity f3 c (h) — 2 /^(Ising), in the regime of 
large \h\, which is expected to vanish in the same limit; 
indeed we have verified that the functional behaviour ex- 
pected from a strong coupling expansion 



f3 c (h) - 2/3 c (Ising) = b ie h + b 2 e 2h + b 



'3C 



3/l 



(17) 



3 This multiplicative factor is caused by a different normalization 
in the Hamiltonians: the Ising one is usually written as a sum of 



fits our data with 61 = 0.0514(1), b 2 = 0.0185(5), 
63 = 0.0150(5) and x 2 /d.o.f. = 2.8. The agreement is 
reasonable taking into account that our data are very ac- 
curate and we truncate the strong coupling series just to 
third order. 

In the opposite limit of small values of h, as shown in 
Fig. [T71 we have been able to fit the /3(h) dependence by 
a third order polynomial in h 



p c {h) = 60 + ^1 h + b 2 h 2 + b 3 h 3 



(18) 



with b = 0.550500(45) 61 = 0.676(44), b 2 = 26(13), 
b 3 = 2100(1400) and x 2 /d.o.f. = 2.5. We notice that 
&o gives an estimate of the critical point position of the 
Potts model without external field compatible with the 
known result (3 c (h = 0) = 0.550565(10) obtained in Q 
and that the slope of j3 c (h) at h — 0~ is different from 
the one observed on the positive h side @ . 

Finally, fitting data for j3 c (h) around h tr ic, we can esti- 
mate also the temperature location of the tricritical point 
and state (/3 tric , h tTic ) = (0.5480(1), -0.00415(3)). 



0.05 - 
0.04 - 



|P 0.03 - 



@ 

u 

co. 



terms ptr. cr,- 



/3(2<5 CT 



1). 



FIG. 16: Plot of p c (h) - 2/3 c (Ising). The line is the result of 
a fit with the function in Eq. (117[) . 



IV. CONCLUSIONS 

We have investigated the critical properties of the 
three-dimensional three-state Potts model as a function 
of a negative magnetic field coupled to one of the three 
spin states. In this case the system possesses a residual 
exact symmetry, which gets spontaneously broken at a 
critical coupling /3 c (h), which approaches twice the criti- 
cal coupling of the 3D Ising model for \h\ — > 00. In par- 
ticular, we have determined the tricritical value /itric at 
which the finite temperature first order transition, taking 
place for null or small values of \h\, turns into a second 
order transition in the universality class of the 3D Ising 
model. 

We have shown that, in proximity of the tricritical 
field /icrit; it is difficult to determine the critical behav- 
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0.55 - 



' 0.549 - 

oo. 



0.548 - 



0.547 I 1 , 1 , 1 , 1 , 1 , L 

-0.005 -0.004 -0.003 -0.002 -0.001 

h 

FIG. 17: Plot of /3 c (h). The line is the result of a fit with the 
function in Eq. [18] The point at h = is the result of || and 
is not included in the fit. The (red in color) spot on the fitted 
line corresponds to our location of the tricritical point. 

ior from the finite size scaling of susceptibilities alone, 
since, at a given distance from h cr n, tricritical scaling 
masks the correct critical indexes up to a given lattice 
size £ max , which is regulated by tricritical crossover ex- 
ponents (I maK oc \h-htric\~ 1 in our case). A better strat- 
egy is to determine parameters which fix the strength of 
the first order region, like the order parameter gap A or 
the latent heat Ae, and to determine h tT i C as the value 
of h at which these parameters extrapolate to zero. We 



also showed that the order parameter cumulant is another 
very useful quantity to look at. In this way we have ob- 
tained the quite accurate estimate htric = —0.00415(3), 
which is in agreement with previous determinations re- 
ported in Refs. [12ll27|. Although in this work the cumu- 
lant method has proved to be the most efficient, which 
of the studied methods is to be preferred to locate a tri- 
critical point should be model dependent. 

Our results may be useful in lattice QCD studies aimed 
at determining the order and universality class of the 
Roberge- Weiss endpoint, and the values of the tricritical 
masses separating the second order from the first order 
regions both for two [1^, |3(| and three [27j degenerate 
flavors. In particular we expect that distinguishing the 
correct critical behavior on feasible lattice sizes will be in- 
creasingly difficult as we approach the tricritical masses, 
the specific heat being anyway more sensitive than the 
order parameter susceptibility. An accurate determina- 
tion of the tricritical masses may be based, for instance, 
on a careful study of the transition strength as a function 
of h in the first order regions. 
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